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We report a quantum Monte Carlo study of the thermodynamic properties of arrays of spin ladders 
with various widths (n), coupled via a weak inter- ladder exchange coupling aj , where J is the intra- 
ladder coupling both along and between the chains. This coupled ladder system serves as a simplified 
model for the magnetism of presumed ordered spin and charge stripes in the two-dimensional Cu02 
planes of hole-doped copper oxides. Our results for n = 3 with weak inter-ladder coupling a — 0.05, 
estimated from the t — t' — t" — J model, show good agreement with the ordering temperature of 
the recently observed spin density wave condensation in La2Cu04+y. We show that there exists a 
quantum critical point at Uc — 0.07 for n — 4, and determine the phase diagram. Our data at this 
quantum critical point agree quantitatively with the universal scaling predicted by the quantum 
nonlinear a model. We also report results on random mixtures of n = 2 and n = 3 ladders, which 
correspond to the doping region near but above 1/8. Our study on the magnetic static structure 
factor reveals a saturation of the incommensurability of the spin correlations around 1/8, while the 
incommensurability of the charge stripes grows linearly with hole concentration. The implications of 
this result for the interpretation of neutron scattering experiments on the dynamic spin fluctuations 
in La2-2!Sra;Cu04 are discussed. 



I. INTRODUCTION 

One of the most surprising early observations in the 
field of high temperature superconductivity is that the 
onset of superconductivity as a function of hole concen- 
tration coincides with a commensurate-incommensurate 
transition in the low-energy dynamical spin fluctuations .El 
These incommensurate spin fluctuations have been stud- 
ied extensively over the past decade.El We note especial 
the recent detailed quantitative study by Yamada et al 
A variety of theoretical explanations of the spin incom- 
mensurability have been offered, varying from n.es.ting of 
the Fermi surface in a nearly free electron modclQ B to mi- 
croscopic phase separation of the holes, ia a doped- Mott 
insulator model for the copper oxides. QlI Specifically, in 
the latter model under certain circumstances the holes 
could organize themselves into stripes with the charged 
stripes acting as antiphase domain walls for the inter- 
vening antiferromagnetic regions. The latter correspond 
to spin ladders which have been mr^ntly studied quite 
extensively in a different context Jl3'EJ None of the ex- 
periments on the spin dynamics has been able to select 



unambiguously between these various theoretical models. 
However, in the past two years, the experimental sit- 
uation has changed significantly. First, rather dramatic, 
elastic magnetic scattering effects have been observed^ 
at low temperatures in samples of Lai.6-xNdo.4Sr2:Cu04 
with a;=0.12, 0.15, and 0.20. These materials are all su- 
perconductors with onset T^s of ~ AK, UK and 15K, 
respectively. In each case, Tranquada and coworkers ob- 
serve elastic incommensurate magnetic peaks with on- 
set temperatjises of ~ 50K, AGK and 15K, respectively. 
Recent workllj has shown that in the a;=0.12 sample 
the correlation length reaches its maximum value below 
^ 30K. The incommensurabilities are essentially identi- 
cal to those of the corresponding dynamical fiuctuations 
in samples of La2-a;Sra;Cu04 with the same x. In their 
study of the spin dynamics in La2-KSr^Cu04, Yamada et 
al. observe that the momentum space width of the low- 
energy spin fluctuations is a minimum for hole concen- 
trations near 1/8. Further, at these hole concentrations 
the spin fluctuations extend down to very low energies 
even in the superconducting state. This is in contrast 
to the situation in samples with hole concentrations at 



and above the optimal doping level, x ~ 0.15, where a 
well-defined spin gap is observed ip-|the superconducting 
state. This has led Suzuki et bLEJ and more recently, 
Kimura et oLEj to search for elastic magnetic scattering 
effects in Lai.88Sro.i2Cu04. Remarkably, they observe 
a spin density wave transition at Tm — SliiT which co- 
incides to within the errors with the onset temperature 
for superconductivity in this sample; the incommensura- 
bility is e = 0.12 which equals the Sr-doping, x. Kimura 
and coworkers also observe magnetic order in a sample of 
Lai.goSrg 10CUO4 (Tj, = ZIK) with incommensurability 
e = 0.105 below Tm o± 15K. _ 

Most recently, Lee and coworkeraij have searched for 
spin density wave order in a sample of La2Cu04 doped 
with oxygen, that is, La2Cu04+i,, which has a supercon- 
ducting Tc (onset) of A2K and is predominantly stage-4. 
They indeed find a transition to long-range incommensu- 
rate magnetic order at Tm — A2K with typical Bardeen- 
Cooper-Schrieffer (BCS) mean field behavior of the or- 
der parameter below Tm- Importantly, they find that 
both the spin ordering direction and the three dimen- 
sional stacking arrangement coincide with those in pure 
La2Cu04. The magnetic structure can be modeled quite 
well with n = 3 stripes of the La2Cu04 in-plane struc- 
ture separated by n = 1 non-magnetic antiphase domain 
walls; this gives an incommensurability e = 0.125, close 
to the measured value. The measured long-range ordered 
moment is ~ 0.15//^ and the inter-layer magnetic corre- 
lation length is ^ 3 Cu02 layers. 

In our view, these observations, including especially 
this recent work on La2Cu04+y, give credence to stripe 
models for the ipicroscopic magnetic structure. The re- 
cent observationEZI of similar incommensurabilities in the 
spin dynamics of underdoped YBa2Cu306.6 strengthens 
this argument. Of course, a great deal of work remains to 
be done to explain all of the experimental observations, 
especially the electronic properties in thet^iormal and su- 
perconducting states. Kivelson et. aZJla consider the 
electronic degrees of freedom within the charge stripes 
and characterize these phases as electronic liquid crys- 
tals. Although the charge and spin degrees of freedom 
of the charged stripes themselves are clearly essential to 
the physics, it is nevertheless important as a first step to 
develop a deeper understanding of the simpler problem 
of the magnetism of idealized insulating stripe arrays as 
a function of stripe width or, equivalently, hole concen- 
tration. Accordingly, we adopt a naive model, in which 
for hole concentrations between ~ 0.05 and Q .15 there is 
approximately one hole per two coppergl3il3 along site- 
centered charge stripes running along the tetragonal Cu- 
0-Cu-O... axes with the holes confined to single chains. 
We assume further that the charge stripes are effectively 
non-magnetic and that the magnetic coupling across the 
stripes is, to first order, determined by the third-nearest- 
neighbor coupling within the Cu02 plane. This is be- 
lieved to be antiferromagnetiCpand to be about 5% of 
the nearest neighbor coupling .E3 We omit entirely any 
effects of the low energy charge and spin excitations of 



the charged stripes themselves. In this greatly simpli- 
fied model, one may use standard quantum Monte Carlo 
(QMC) techniques to determine the magnetic properties 
of the stripe arrays. 

Our understanding of two-dimensional quantum mag- 
netism has grown enormously since the discovery of 
high temperature superconductivity. Theoretical ef- 
forts to understand the underlying antiferromagnetism 
in the parent compounds as well as experimental ef- 
forts to synthesize and study new cuprate materials 
have both contributed significantly. In particular, the 
two-dimensional quantum Heisenberg antifcrromagnet 
(2DQHA) has been studied extensively, since the pioneer- 
ing work by Chakravarty, Halperin, and Nelson.Eil They 
have mapped the long- wavelength, low-temperature be- 
havior of the 2DQIIA to a quantum nonlinear a model 
(QNLcrM) in (2-1-1) dimensions and have obtained a 
phase diagram with three regimes; quantum disordered 
(QD), quantum critical (QC), and renormalized classical 
(RC). They have argued that the 2DQIIA on the square 
lattice for S > 1/2 with nearest neighbor interactions 
is in the RC regime, and consequently has a correlation 
length diverging exponentially in 1/T as temperature is 
lowered to T = 0, implying the existence of a long-range 
ordered ground state. The temperature dependence of 
the correlation length in the RC regime, which has been 
solved exactly to-three loop order by Hasenfratz and Nie- 
dermayer (IIN)J23 agrees quantitatively with the results 
of neutron scattering experiments on both Sr2Cu02Cl2E^ 
and La2Cu04.ca 

One-dimensional (ID) quantum magnetism has also 
drawn much interest in the last decade, mainly due to 
a conjecture made by Haldane regarding the different 
ground state properties between spin chains with inte- 
ger and half-integpe spin quantum number .E3 In his sem- 
inal work in 1983J23 Haldane mapped the ID Heisenberg 
model onto the QNLcrM and showed that the half-integer 
spin chain has an additional topological term, which gen- 
erates gapless low-energy excitations together with al- 
gebraically decaying spin correlation functions. Integer 
spin chains, however, can be described by the standard 
QNLcrM in (H-l) dimensions, in which the ground state 
is disordered due to quantum fluctuations and the exci- 
tation spectrum acquires a gap - the so-called Haldane 
gap. This conjecture has been subsequeatly confirmed 
both numericallyEII and experimentallyBa Spiui-Jadders 
have been studied mostly in the same context ;t2l when 
even numbers of S' = 1/2 chains are coupled to form a 
ladder (even- width) , they show the same universal be- 
havior as integer spin chains, while odd-width S ~ 1/2 
ladders behave essentially like a single 5 = 1/2 chain at 
low temperatures and long wavelengths. 

In this paper we report a detailed QMC study of ladder 
arrays as a function of the width of the spin ladders (n) , 
the strength of the coupling across the non-magnetic line 
(J' — aJ), and temperature (T). We have carried out 
such calculations for n = 4 ladder arrays, n = 3 ladder 
arrays, and for random mixtures of n = 3 and n = 2 lad- 



ders. These various arrays are illustrated in Fig. |l|. In the 
limit of a — > 0, the model corresponds to a set of isolated 
ladders, and in the limit a — > 1 we recover the isotropic 
S — 1/2 square lattice QHA. One can easily see that a 
quantum critical point as a function of a should exist for 
non-zero a for coupled even-widtii. ladders. Indeed, the 
recent study by Tworzydlo et alBB shows that ac ~ 0.3 
for coupled n = 2 ladders. We show that arrays of n = 4 
ladders exhibit a quantum critical point around a ~ 0.07. 
We discuss our QMC results in the context of available 
experimental information; we also suggest future mea- 
surements which should determine whether or not these 
calculations are in fact relevant to the real monolayer 
Cu02 superconductors. 

The format of this paper is as follows: In Sec. H we 
discuss our QMC techniques. Section III contains our 
results for weakly coupled three-leg ladders (n — 3). Re- 
sults for weakly coupled four-leg ladders (n = 4) includ- 
ing especially the quantum critical behavior are given in 
Sec. [V, In Sec. 0, we discuss the results for random 
mixtures of two-leg and three-leg ladders. We have also 
studied strongly coupled n = 3 ladders and present the 
Finally, a discussion, summary, and 



VI 



results in Sec. 

conclusions are given in Sec. VII 




FIG. 1. Schematic diagram of the model, (a) Weakly cou- 
pled three leg ladder array with a periodicity of 8a. J' is the 
inter-ladder coupling aj. (b) Weakly coupled four leg ladder 
array, (c) Random mixture of two and three leg ladders. 



II. QUANTUM MONTE CARLO 

We have carried out quantum Monte Carlo simulatioBS, 
on large lattices utilizing the loop cluster algorithm.E2l 



The sapap algorithrrip^reviously employed to study spin 
ladderalj and chainsEil is used with minor modifications. 
We obtain the temperature dependence of the uniform 
susceptibility, Xu(T), the spin-spin correlation length in 
the x-direction, ^x{T), and the y-direction, ^y{T), the 
staggered susceptibility, Xs{T), and the static structure 
factor at the antiferromagnetic wave vector Q = (7r,7r), 
Sq{T). The Hamiltonian for our model is based on the 
square lattice with a coupling J, but every n-th bond 
in the x-direction is replaced by a bond with exchange 
coupling a J: 



n = J 



/ . ^hj ■ Sj.j+i + } ^ 



ii+i.j 



*j 



i^ln,j 



i—ln,j 



5j+l,j 



(1) 



where i and j run over the x and y directions, respec- 
tively, n is the width of the ladder and I counts the 
ladders. We use units in which h = Ub = QIJ-b = 1- We 
also set the lattice constant a — \. Periodic boundary 
conditions are used in both directions. 

One should note that, although this Hamiltonian cap- 
tures the essential magnetic interactions of coupled lad- 
ders, if Eq. (P) is to describe the stripe arrays shown in 
Fig. P, the lattice constant in the x-direction is -^^ in- 
stead of one. This affects only two quantities in our com- 
putation; ^x{T) and Sq{T). The correction for ^x{T) is 
trivial, since £,x{T) is measured in units of the lattice con- 
stant. Unless noted otherwise, this ^fr correction is in- 
corporated in the ^x data presented in this paper. Sq{T), 
however, cannot be converted in a simple manner for the 
stripe model. Therefore, in the following sections, Sq{T) 
should be understood as the structure factor at (tt, tt) 
for the model Hamiltonian, Eq. (|l|), rather than for the 
stripe model as shown in Fig. nl The incommensurate 
structure factor is discussed only in Sec. 0, where it is 
separately computed as a function of Qx, S{qx, tt). 

In order to deduce the spin-spin correlation lengths, 
^xiT) and ^y{T), the instantaneous spin-spin correlation 
function, C{x,y), is computed and fitted to the asymp- 
totic form 



C(x,0) - a;-^e-^/«- + {Lx - x) 



-Ap-(ix-a^)/Cx 



C{0,y) 



y~^e~y/^y 



(Ly 



(2a) 
y)-Ag-(L„-a)/««^ (2b) 



which is a symmetrized 2D Ornstein-Zernike form (A = 
0.5). Only data with x ^ 3^2^ and y ^ 3^j^ are included in 
the fits to ensure that the asymptotic behavior is probed. 
Technically, due to the non-uniformity of the exchange 
couphngs in the x-direction, C{x,Q) can be better de- 
scribed by the function ^ x^^'^e"^'^"^ sin(a;/n) rather 
than Eq. (Ea); however, the values for ^x obtained with 
both methods agree within the error bars. 

The lengths and Trotter numbers of the simulated lat- 
tices are chosen so as to minimize any finite-size and 



lattice-spacing effects. The linear sizes of the lattice, L^ 
and Ly, are kept at least 10 times larger than the respec- 
tive correlation lengths. Spin states are updated about 
2 X 10"* times to reach equilibrium and then measured 
~ 5 X 10* times. For the random mixtures of n = 2 and 
71 = 3 ladders, we generate each configuration in the fol- 
lowing way: Beginning from the first column {i — 1), we 
assign a spin or hole for each column; the number of spin 
columns between the hole columns represents the width 
of the spin stripe. In order to make sure that the width 
is either two or three, we put restrictions such that a spin 
column must be assigned after each hole column, while 
a hole column must be assigned after three consecutive 
spin columns. The width of each ladder is determined by 
generating a random number r (0 < r < 1) and compar- 
ing r with p, where p is the ratio of the three leg ladder 
in the mixture. After two spin columns are assigned, the 
spin column is assigned only if r < p, otherwise the hole 
column is assigned and the width becomes two. For ex- 
ample, for an equal mixture of two and three leg ladders, 
p = 0.5. For each T and p, typically 10 to 20 different 
configurations are generated and averaged over. Averag- 
ing over more configurations does not alter our results. 



III. WEAKLY COUPLED THREE LEG LADDERS 

We present our QMC data for the n — i case. Fig. 
Il|(a), inj-this section. From detailed studies of isolated 
ladders ,liij it is known that at low temperatures and long 
length scales three-leg ladders exhibit the same behavior 
as a single chain [S — 1/2). It is well known that any 
non-zero inter-ladder coupling eventually enhances the 
correlations across the non-magnetic stripe and drives 
the system towards the 2DQIIA-QNLctM RC fixed point. 
Thus, we expect to see qualitatively different behavior for 
an array of weakly coupled n = 3 ladders from that of 
isolated ladders. The correlation lengths for arrays of 
n = i ladders are shown in Fig. as a function of inverse 
temperature for various inter-ladder couplings {a). Note 
that a = corresponds to an isolated three leg ladder; 
these data are taken from Ref. O. Not surprisingly, one 
can see clearly ^ exp(l/r) behaviors except for a = 0, 
as discussed above. In the inset of Fig. ||, we also show 
the ratio £,y/£,x- At low temperatures, ^y/^x approaches 
the mean field prediction, y/a, shown here as a solid line 
for each a. 

Our correlation length data are fitted-io the crossover 
form given by Castro Neto and Hone,E3 which interpo- 
lates between the HN result at low temperatures and 
^ ~ T~^ at high temperatures: 



^ = Aexp(— ^) 



T 'I 



1_T_ 
2 2irps 



(3) 



where Ps is the spin stiffness. Our fitting results are pre- 
sented in Table H. We note first that the values deduced 
for 2'Kps separately from ^j, and ^^ are the same within 



the error bars, which means that there is only one tem- 
perature scale for the low-temperature behavior of this 
model, despite the fact that the correlation lengths them- 
selves are highly anisotropic. 
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FIG. 2. Correlation lengths for arrays of n = 3 ladders 
with different inter-ladder coupling a plotted as a function of 
inverse temperature. Correlation lengths in the x-direction 
are shown as solid symbols, while those in the y-direction are 
shown as empty symbols. For a = 0.05, the fitting results 
from Eq. (H) are also plotted. Inset: The ratio of ^y/^x- Each 
solid line corresponds to the mean field value, ^/a. The ratio 
is given without the n/{n + 1) correction discussed in Sec. Ill 



TABLE I. Fitting parameters for n=3 data to Eq. (H) with- 
out the n/{n + 1) correction. 



2-Kps/J 
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iv 
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0.05 0.10(1) 0.11(1) 

0.1 0.17(2) 0.19(1) 

0.3 0.55(2) 0.58(1) 

l.O'' 1.1310(3) 



0.26(2) 1.14(5) 

0.52(5) 1.4(1) 

0.83(5) 1.22(5) 

0.4978(8) 
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FIG. 3. Uniform susceptibility, staggered susceptibility, 
and inverse static structure factor peak intensity for arrays 
of n = 3 ladders with different inter-ladder couplings a. 

Since the 2D Heisenberg model cannot have long- 
range order at any non-zero temperature, either a small 
anisotropy in the spin Hamiltonian or a non-zero 3D 
coupling is necessary to explain the Neel transition in 
real materials. In pure La2Cu04, one can estimate the 
Neel transition temperature reliably as the temperature 
at which the correlation length squared becoipes of or- 
der of the inverse of the effective anisotropyJ23 which, 
in the case of La2Cu04, is about 7e// ~ 10^^. We 
can apply the same heuristic formula to determine the 
spin density wave ordering temperature, T™, for our pre- 
sumed stripe array model for doped La2Cu04. For ex- 
ample, La2Cu04+j^ has the same spin ordering direction 
and stacking scheme as in pure La2Cu04 except that the 
inter-layer order extends over only about three planes;E£l 
thus, we take the same effective anisotropy and find the 
temperature satisfying £,x{Tm)iy{Tjn) ~ 77//- In our 
simplified model the inter-ladder coupling is determined 
by the third-nearest-neighbor exchange coupling in the 
Cu02 plane; from curxent values for t, t' ,t" and J in the 
extended t — J modeltj we estimate a « {t" /tY ~ 0.05. 
By extrapolating our correlation length data to lower 
temperature for a — 0.05, we obtain Tm ~ 0.029J, or 
Tm « 44iir for J = 1500iir, which agrees remarkably well 
with the experimental result. This precise agreement is 
a coincidence since a = 0.05 is only a crude estimate and 



we have assumed perfectly ordered charge stripes. We 
note that a = 0.1 gives Tm ~ 7hK. We should also note 
that these estimates of T,„'s are not very sensitive to the 
explicit value of 7e//- Even if 7e// changes by a fac- 
tor of two, Tm remains within 10% of the value obtained 
here. However, Tm depends sensitively on the choice of 
the inter-ladder exchange coupling aJ . The low values 
for Tm found experimentally clearly constrain the inter- 
ladder coupling to rather small values (or, as discussed in 
Sec. ^, very large values), at least within the simplified 
model considered here. 

In Fig. |3[ Monte Carlo results for x„, Xs, and the in- 
verse of Sq are plotted as a function of T, for different 
a's. The reason Sq^ is plotted is to emphasize the effect 
of the inter- ladder coupling. One of the most significant 
results of Ref. ^ is that the divergence of Sq is only log- 
arithmic in ID; thus Sq does not extrapolate to zero at 
zero temperature. As shown in Fig. [^(c), coupled n = 3 
ladders behave completely differently from the isolated 
n = 3 ladder; namely, the Sq of coupled ladders exhibits 
a clear crossover from the weak divergence in T^^ of an 
isolated ladder to the strong divergence in T~^ of a 2D 
spin system. 
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FIG. 4. Phase diagram for weakly coupled four leg ladders. 
The quantum critical point corresponds to a = 0.07(1), as 
shown in Fig. M. The shaded line represents a crossover region 
rather than a true phase boundary. 



IV. WEAKLY COUPLED FOUR LEG LADDERS 



Results for n = 4, corresponding to Fig. |l|(b), are pre- 
sented in this section. Unlike the n = 3 case, an isolated 
four-leg ladder has a Haldane gap (A k. 0.16 J). There- 
fore, if aJ <C A, the ground state remains disordered, 
and there must then be a quantum phase transition as a 
function of a at zero temperature at a critical value ac to 



an ordered ground state. In QNLcrM language, there is a 
quantum critical point dividing the QD ground state and 
the RC ground state at a^- We show the phase diagram 
of arrays of n = 4 ladders in Fig. |j. In this figure, we 
take lixps as the temperature scale where crossover from 
the QC to the RC regime occurs, and A as the crossover 
temperature from the QC to the QD regime. I'Kps is 
determined in the manner explained in Sec. Ill; A can 
be obtained by fitting to Xu{T) ~ exp(— A/T) at low 
temperatures (T ^ A). One should note that the thick 
shaded lines in Fig. represent crossover lines rather than 
true phase boundaries. The pihase diagram shown here is 
basically that of the QNLcrM,EJ with a playing the role of 
the coupling constant g. As shown in Fig. ^, a quantum 
critical point occurs at a^ — 0.07(1). 
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FIG. 5. Uniform susceptibility, staggered susceptibility, 
and inverse static structure factor peak intensity for arrays 
of n = 4 ladders with different inter-ladder couplings a. 

Monte Carlo results for Xu, Xs: and S'q^ are plotted as 
a function of T, for different a's, in Fig. ^ The corre- 
lation length obtained from the simulation is plotted in 
Fig. H as a function of T^^. Figures ^ and g are plotted 
in such a way as to contrast their behaviors with those of 
arrays of n = 3 ladders. The difference between the be- 
haviors in the QD and the RC regime is evident in these 
figures; namely, ^, Xsi and Sq all diverge exponentially 
at low temperatures in the RC regime (e.g., a = 0.3), 
while they all saturate at finite values as T ^ in the 



QD regime (e.g., a = 0.05). 
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FIG. 6. Correlation lengths for arrays of n = 4 ladders 
with different inter-ladder couplings a plotted as a function 
of inverse temperature. Correlation lengths in the a;-direction 
are shown as solid symbols, while those in the y-direction are 
shown as empty symbols. 

To illustrate the quantum critical behavior more dra- 
matically, we plot in Fig. R(a), the dimensionless ratio 
Sq/Txs- At very high temperatures, this ratio shows 
classical behavior, Sq/Txs = 1, while the same behav- 
ior shows up again at very low temperatures for a> cuc, 
that is, in the RC regime. This is as expected, since 
RC behavior is closely similar to classical behavior, if the 
spin wave velocity and spin stiffness renormalizations are 
accounted for. In the QD regime, that is, a < otc-, Sq 
and Xs are constant; therefore Sq/Txs should be linear 
in T~^ . According to thfLxiuantum critical scaling pre- 
dicted for the QNLcrMpl'L^ this ratio should show uni- 
versal behavior in the QC regime with the specific value: 
Sq/Txs = 1.10(2). As may be seen in Fig. 0(b), the 
a = 0.07 data for Sq/Txs indeed are constant ~ 1.12 
at low temperatures, in quantitative agreement with the 
QC theory. This strongly supports the above claim that 
ftc — 0.07 for the n = 4 stripe array system is a quantum 
critical point. Another quantity plotted in Fig. ^(b) is 
S,yT. Although ^xT is not shown here, it shows essen- 
tially the same temperature dependence. Again at very 
high temperatures, this quantity shows classical behav- 
ior, ^j,T = J, while the low-temperature behavior can 



distinguish the QC regime froni|tiie RC and QD regimes. 
In the QC regime, the QNLcrMld predicts £T = c/1.04, 
where c is the spin wave velocity. As may be seen in Fig. 
0(b), S,yT is indeed a constant at low temperatures (high 
1/T) for a = 0.07; the explicit value ^yT ~ 1.4J corre- 
; 1.46 J, close to, but somewhat less than. 



sponds to c 

the value c — 1.657 J for the square lattice.E^ 
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FIG. 7. (a) The ratio Sq/Txs as a lunction of inverse tem- 
perature as described in the text, (b) The correlation length 
in the y-direction of arrays of n = 4 ladders multiplied by tem- 
perature is plotted to contrast the different low-temperature 
behaviors in the QD and RC phases. The solid lines in both 
(a) and (b) are the predictions for quantum critical scaling. 



V. RANDOM MIXTURE OF WEAKLY 
COUPLED TWO AND THREE LEG LADDERS 

In Fig. Ma), we show representative static structure 
factor data at low temperatures for random mixtures of 
n = 2 and n = 3 weakly coupled ladders, obtained as de- 
scribed in Sec. O. Here, p is the fraction of n = 3 ladders 
in the mixture; therefore, p — 1.0 corresponds to the pure 
n — 3 case. We choose a — 0.05 for all of the simulations 
discussed in this section. We show data at low temper- 
atures (T ~ 0.05 J), since the magnetic structure factor 
develops well defined peaks only at these low tempera- 
tures. At high temperatures, due to the short correlation 
lengths, the structure factor exhibits only broad peaks, 
making it difficult to extract meaningful values for the in- 
commensurability. Nevertheless, we can fit the data with 



two identical Lorentzians, split symmetrically about the 
antiferromagnetic wave vector Q = (7r,7r), together with 
a broad temperature-independent background term cen- 
tered around {tt,tt): 



1 



1 



'^^'^"'^^" 4 U + ql/^^l ' l + g2/^2 
B 



1 + (fe - 7r)V 



2/^2- 



(4) 



Here q± = Qx — ('"'ie) and the last term is a temperature- 
independent background. The width of this background 
term corresponds to ~ 1 lattice constant. There are only 
two temperature dependent fitting parameters, e and Kx, 
since the peak intensity Sq/A is calculated separately in 
our QMC study. The peak position corresponds to the 
incommensurability e, and the peak width corresponds 
to the inverse correlation length Kx ■ These data are pre- 
sented in Fig. ^. In Fig. ^(b), we also show Ky = 
S,y is computed directly in the way explained in Sec. IH 
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FIG. 8. (a) Static Structure factor near the antiferromag- 
netic wave vector for the n — 2 and n — 3 mixed spin ladders. 
Representative data at T = 0.05 for p — 0.5 and p — 0.25, 
and at T = 0.08 for p = 1.0 are shown; a = 0.05 for all 
data. The solid lines are the results of fits to two Lorentzians 
[Eq. (0)]. (b) Corresponding charge structure factor near the 
nuclear zone center, calculated as described in the text. We 
show p = 0.95 data instead of p = 1.0 for graphical purpose. 



The most prominent feature in Fig. |8|(a) is that the in- 
commensurabiHty e of the magnetic structure factor does 
not change significantly as p is varied. This is verified 
quantitatively in Fig. y(a), where it may be seen that the 
values for e for both p — 0.5 and p — 0.25 cluster around 
1/8, the exact value for p — 1.0. Note that the high- 
temperature data have large error bars, since the peaks 
are very broad. This behavior of the magnetic structure 
factor may be contrasted with that of the charge struc- 
ture factor. We show in Fig. ||(b) the charge structure 
factor for the charge stripes in the random mixtures. The 
calculation is made for the simplest charge distribution: 
we assume that the scattering is unity from the charge 
stripes (anti-phase domain walls) and zero from the spin 
ladders. Although greatly oversimplified, as an illustra- 
tion this calculation nevertheless provides useful intuitive 
guidance. The most important qualitative result is that 
the incommensurability of the charge stripes is not equal 
to twice the magnetic incommensurability in the random 
mixtures. For example, the p = 0.25 case has a charge 
incommensurability of 0.31, which is quite different from 
2e — 0.25(2), where e is the magnetic incommensurabil- 
ity. 
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FIG. 9. (a) The incommensurability obtained from the fit- 
ting static structure factor to two Lorentzians together with 
a fixed background term [Eq. (W)] . Even if there is a substan- 
tial fraction of n = 2 ladders, the incommensurability remains 
fixed near 1/8. (b) Inverse correlation length as a function of 
T. Kx is obtained from the fit. K,y is determined from the 
correlation function Eq. (Ehl) as described in Sec. |ll 

This saturation of e of the magnetic structure factor 



around 1/8 can be understood in the following way: At 
low temperatures (T <C A, A « 0.41J) the spins on 
the two-leg ladder will form a spin singlet ground state 
with a large spin gap in the excitation spectrum. These 
spin singlets have effective spin zero, and therefore do 
not contribute in first order to the magnetic structure. 
Therefore, the magnetic structure factor originates pre- 
dominantly from the n = 3 components of the mixture. 
In their systematic study of the dynamic incommensu- 
rate spin Jluctuations in La2-3;Srj;Cu04, Yamada and 
coworkersn showed that the incommensurability is linear 
in hole concentration x with e ~ a; for 0.06 ^ a; ^ 0.12 
and saturates around 1/8 on further doping. In our sim- 
plified model of stripes, the range 0.12 ^ a; < 0.16 cor- 
responds to mixtures of n = 2 and n = 3 ladders. Al- 
though we have calculated the static structure factor, we 
believe that the saturation in the dynamic fluctuation in- 
commensurability in La2_a;Sra;Cu04 may have the same 
origin. We also have computed Xmj Xs; and Sq , which 
are plotted in Fig. O as functions of T, for different p's, 
all with a = 0.05. Note the differing low-temperature 
behaviors of the p = arrays from those with p ^ 0; it is 
evident that the n = 3 physics indeed dominates in the 
magnetism of the random mixture. 
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FIG. 10. Uniform susceptibility, staggered susceptibility, 
and inverse static structure factor peak intensity for arrays 
of n = 2 and n = 3 mixed ladders with different p. 

Finally, we should note that in our model, in order to 



explain the fact that the spin incommensurability sat- 
urates near 1/8 even for very high dopings, we must 
hypothesize that for doping above the optimal value 
X ~ 0.15, the charge per stripe increases progressively 
with increasing x above ~ 0.15. That is, we assume that 
as the doping increases beyond x = 0.15 charge stripes 
cannot become closer due to the Coulomb repulsion; in- 
stead, additional charges go into already existing charge 
stripes. 



VI. STRONGLY COUPLED THREE LEG 
LADDERS 

The previous calculations have all been done in the 
weak inter-ladder coupling limit; indeed, in order to ob- 
tain reasonable magnetic ordering temperatures in our 
model, the reduced inter-ladder coupling must be of otdjer 
a ~ 0.05. However, as pointed out to us by Kivelson,EZI a 
similar situation also should obtain in the limit of large 
a. In this case, physically one would imagine the car- 
riers along the charge stripes mediating a large effective 
antiferromagnetic exchange between the bounding spin 
chains; these neighboring chains would then form an ar- 
ray of two-leg ladders coupled to single chains [Fig. |l|(a)]. 

Accordingly, we have carried out a limited number of 
calculations on arrays of n = 3 ladders coupled with 
strong exchange couplings; that is, in the a > 1 limit. 
Specifically, we have repeated our simulations of Sec. Ill 
with a = 4.0 and a = 10.0. The correlation length data 
so-obtained are shown in Fig. 0(a). To facilitate a com- 
parison, we have also plotted the a = 0.05 data from 



Fig. 2. At low temperatures, one can see the RC behav- 
ior clearly, similar to that of the weak coupling data. The 
solid lines are the results of fits to Eq. (||). It is evident, 
however, that the uniform susceptibility data for strongly 
coupled n = 3 ladders shown in Fig. |l l|(b) are quite dif- 
ferent from those of weakly coupled {a = 0.05) n = 3 
ladders. In this limit, since J' is the primary coupling, 
as discussed above, two-leg ladders are effectively formed 
across the charge stripes; therefore, we should consider 
this system as a mixture of 5 = 1/2 chains and n — 2 
ladders rather than as an array of n = 3 ladders. Specif- 
ically, S = 1/2 chains and n = 2 ladders alternate in the 
x-direction. The spin gap of these n — 2 ladders is large 
(of order 0.5 J') and the contribution from these to Xu is 
essentially zero. Thus one might surmise that Xu would 
be entirely due to the S — 1/2 chains. That this idea is 
correct may be seen from Fig. pT|(b) , where we have plot- 
ted the uniform susceptibility per spin of the S = 1/2 
chain from Ref. 31 divided by three. Both the a = 4.0 
and the a = 10.0 data agree with the spin chain results 
for T ^ 0.4J. At low temperatures, the correlations be- 
tween the chain and the ladder become important and 
we observe behavior similar to that of mixtures of lad- 
ders [Fig. |l^(a)]. 

The correlation length for a — 4.0 grows as a function 
of T^^ much faster than that of a = 0.05, such that the 



prediction for Tm, in this case is Tm « lAAK for J — 
1500iir, which is more than a factor of three larger than 
the experimental result. In order to obtain reasonable 
values for Tm in the large a limit, one needs to have 
a close to 10 as shown in Fig. O. In such a case, for 
a — 10.0, we obtain T,n ~ 55K. Such large values for 
J' w 1.3eV seem to us to be unlikely. For this model 
(J' ^ J) to be relevant one would need to invoke some 
other mechanism such as extensive structural disorder or 
effects from charge and spin fluctuations of the charged 
stripes themselves to reduce Tm down to the observed 
values of 30K to iOK. 
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FIG. 11. (a) Correlation lengths for arrays of n = 3 lad- 
ders with strong inter-ladder coupling (a > 1) plotted as a 
function of inverse temperature. The correlation lengths in 
the K-direction are shown as solid symbols, while those in the 
j/-direction are shown as empty symbols, (b) Uniform suscep- 
tibility per spin as a function of T. 



VII. DISCUSSION 

A few comments on the applicability of our model are 
in order. First and foremost, our model is clearly an 
extremely simplified version of any actual charge and 
spin ordering in the doped copper oxide planes. We as- 
sume that (1) the charge stripes are perfectly ordered 



with the charges confined to Cu-O-Cu-0 chains and (2) 
the charge and spin degrees of freedom along the stripes 
can be ignored. We know that dynamic transverse fluc- 
tuations of the charge stripes could be very large, espe- 
cially for incommensurate charge stripe spacingsJ13 How- 
ever, this effect may be less significant if the charge or- 
dering occurs at a much higher temperature than the 
spin ordering, so that the charge fluctuations are signifi- 
cantly reduced near the spin ordering temperature. Sev- 
eral neutron scattering experiments on La2-a;Sra;Ni04, 
and Lai.6-a:Ndo.4Sra;Cu04 indeed report a charge or- 
dering temperature that is higher than the spin order- 
ing temperature.LS The charge degrees of freedom in the 
metallic stripes, which we have ignored, clearly are essep. 
tial for the transport and superconducting properties,E3 
and we have nothing to contribute on this aspect of the 
problem. Second, microscopic phase separation of holes 
and spins can occur in other geometries; for example, it 
is possible to have grid-like spin-rich regions separated 
by hole-rich domain walls as well as diagonal stripes. In 
both cases, the incommensurate spin density wave peaks 
would be rotated by 45° with respect to those for the 
stripe model discussed here. In a recent experiment on an 
X = 0.05 sample]^ which notably is an insulator rather 
than a superconductor at low temperatures, 45° rotated 
magnetic peaks are indeed observed; in that case it is 
believed that this rotation is due to the establishment of 
diagonal stripes. 

Despite the success of the current model in explain- 
ing a number of experimental observations for incom- 
mensurabilities near 1/8, many points remain to be un- 
derstood. As is obvious in Figs. ||, Q and ^ ^^ and 
^y are very anisotropic, which would, in turn, predict 
anisotropic widths in the quasielastic and dynamic neu- 
tron scattering measurements. Currently available exper- 
imental data are not complete enough tp-test this predic- 
tion. In recent work, Tranquada et a/Ji3 calculated the 
elastic structure factor oi n — 4: and n = 3 mixtures 
for sinusoidally varying charge and spin density waves; 
they hypothesized that the stripe disorder is caused by 
pinning by the random Sr dopants. They also argued 
that these defects would explain the absence of any peak 
width anisotropy in their data. 

An obvious major deficiency of the current model is 
that it does not naturally explain the special role which 
hole dopings x and concomitantly, incommensurabilities 
e near 1/8 play in the La2Cu04-based superconductors. 
It is evident that T,„ will be a maximum locally around 
x ~ e ~ 1/8, since for small a, admixtures of n = 4 
{x < 1/8) or n = 2 (x > 1/8) ladders will decrease 
the magnetic correlation lengths and hence will lower 
Tm- However, there is no obvious reason why other odd- 
width ladder arrays would not also be favorable. We 
note that for x. < 0.05 = 1/20 the spin correlations are 
commensurateErB so that the uniform stripe model does 
not apply. Thus, the other relevant odd- width ladder ar- 
rays are n = 7, corresponding to a; = e = 0.0625, and 
n = 5, corresponding to a; = e = 0.083. Wc have carried 



out some limited calculations for n — 5|-and, as one would 
intuit from the results of Greven et aZ.EJ for isolated lad- 
ders, at a given temperature the correlation lengths for 
arrays of n = 5 ladders are somewhat larger than those 
for n — 3 arrays. For perfectly ordered charge stripes, Tm 
should be correspondingly higher. In fact, the opposite 
seems to be true, although more experimental data are 
required to document this completely. At the minimum, 
from Ihe inelastic neutron scattering studies of Yamada 
et aL,Q we know that the low-energy dynamical coherence 
length at low temperatures for a;~e~l/12is shorter 
than that for x ~ e ~ 1/8. 

Clearly, therefore, if this stripe model is to describe 
the real La2Cu04-based superconductors, then some fur- 
ther physics is required. At least three different factors 
could act to reduce the charge and concomitantly the 
spin correlation lengths at lower doping. First, as evi- 
denced by the commensurate spin structure for x ^ 0.05, 
pinning of the stripes becomes much more important at 
lower dopings. Second, as the charge stripe separation 
grows, the stripe-stripe interaction energy will decrease 
and, correspondingly, stripe positional fluctuations will 
grow. These would in turn inhibit the development of 
the spin correlations. Third, as discussed above, it is 
believed that at a; = 0.05 the stripes switch from being 
approximately along the Cup£)-Cu-0 direction to being 
along the diagonal direction,L£l that is from (1 0) to (1 1). 
Presumably, as the doping decreases towards 0.05, fluc- 
tuations into the diagonal stripe phase could occur and 
thus would shorten the spin correlation length. These 
ideas are, of course, purely heuristic. Some real theo- 
retical calculations plus further experiments are clearly 
required to put the model on a firmer basis. 

We have assumed that the charge stripes are site- 
centered rather than bond-centered.c2l The basic reason 
for this is that, at least within our simplified model, we 
are unable to match at all the observed experimental 
trends with a bond-centered model. This may be eas- 
ily seen as follows. First, if the spins across the charged 
stripe bonds are coupled antiferromagnetically, then the 
basic commensurate antiferromagnetic situation is unal- 
tered. Second, if, as originally argued by Aharony et 
al.^3 holes on the oxygens induce a strong ferromagnetic 
coupling between the two neighboring copper spins, then 
one has effective S = 1 chains. Thus, for the specific 
case of e = 1/8, the system breaks up into alternating 
spin-1 chains and two- leg ladders. Both the spin-1 chain 
(A ss 0.4 J) and the two-leg ladder (A « 0.5 J) have large 
spin gaps and short correlation lengths, so that this sys- 
tem only orders at extremely low temperatures, if at all. 
On the other hand, for x = 1/10 or x = 1/6, one has 
single chains or three-leg ladders respectively in between 
the spin-1 chains. Hence the correlation length along the 
chain diverges at low temperature, making a spin density 
wave ordering possible. The difficulty for a bond-centered 
model then is that one has facile ordering at x — 1/10 
and X = 1/6, but not at 1/8, which is the opposite of 
the general trend observed experimentally. Therefore, 
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we are required within the context of our model to adopt 
site-centered stripes in order to obtain sensible behavior 
around x — 1/8 doping regime. It would clearly be in- 
teresting to see if this same difficulty occurs iof the t~ J 
model as discussed by White and Scalapino.E3 

Manifestly, more experiments are needed to address 
many unresolved features of the stripe model. First, one 
needs to understand in detail the doping dependence of 
the spin density wave ordering in the La2-xSr2;Cu04 and 
related systems. Second, better characterization of the 
charge stripes themselves, including both the incommen- 
surabilities and the correlation lengths is crucial. Third, 
unambiguous determination of the charge ordering tem- 
peratures is essential, since it is still not certain as to 
whether charge ordering drives the spin ordering, or vice 
versa.c3 Since neutrons scatter from the small nuclear 
displacements induced by the modulated charge density, 
the charge ordering peaks observed via neutron scatter- 
ing are extremely weak. Electron diffraction is, of course, 
especially sensitive to the charge modulation, but it is 
only useful for surfaces or very thin films which may dif- 
fer from their bulk counterparts. X-ray scattering seems 
to be a logical choice to study the charge modulation in 
these materials. In their study of Lai.48Ndo.4Sro.i2Cu04, 
Zimmermann et al£3 used high-energy x-ray (100 keV) 
to investigate charge stripes, |Ganfirming results from neu- 



tron scattering experiments. O They were able to deter- 
mine peak widths at various temperatures with some- 
what better precision than that of neutron scattering; 
however, a much higher resolution x-ray study is nec- 
essary to compare quantitatively the correlation length 
data with various theoretical predictions. Generally, it is 
important to determine the nature of the charge-charge 
stripe correlations in both the normal and superconduct- 
ing states. 

In summary, we have studied the magnetism aris- 
ing from the charge and spin stripe order in monolayer 
cuprate superconductors. Quantum Monte Carlo simula- 
tions have been carried out on a simplified model, which 
is essentially weakly coupled insulating spin ladders. Our 
calculations are consistent with recent results of spin den- 
sity wave ordering and dynamical spin fluctuations in 
La2_a;Sra;Cu04 and La2Cu04+y for doping and incom- 
mensurabilities near 1/8. We show that the periodicity 
of the incommensurate spin order is not necessarily twice 
that of the charge order. The behavior at lower doping 
remains problematic. We have also studied in detail the 
quantum critical behavior in coupled four-leg ladders. 
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